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FORTRAN PROGRAM FOR PREDICTING OFF-DESIGN PERFORMANCE 
OF CENTRIFUGAL COMPRESSORS 
by Michael R. Galvas 


Lewis Research Center and 
U. S. Army Air Mobility R&D Laboratory 

SUMMARY 

A FORTRAN program for calculating the off -design performance of centrifugal com- 
pressors with channel diffusers is presented. Use of the program requires complete 
knowledge of the overall impeller and diffuser geometries. Individual losses are com- 
puted using analytical equations and empirical correlations which relate loss levels to 
velocity diagram characteristics and overall geometry. On a given speed line com- 
pressor performance is calculated for a range of inlet velocity levels. At flow rates 
between surge and choke, individual efficiency decrements, compressor overall effi- 
ciency, and compressor total pressure ratio are tabulated. An example case of per- 
formance comparison with a compressor built by a commercial engine manufacturer is 
presented to demonstrate the correlation with limited experimental data. 


INTRODUCTION 

Centrifugal compressors are frequently used in refrigeration cycles, refining sys- 
tems, aircraft auxiliary systems, turboshaft engines, and other systems where light- 
weight compact compression is required. The off -design performance characteristics 
of centrifugal compressors are of interest because of the large effects that compressor 
component performance has on overall cycle performance and because the compressor 
is required to operate at off -design conditions much of the time. For turboshaft engines 
assuming a drive turbine inlet temperature of 1422 K and efficiency of 90 percent, a 
2.6-percent loss in compressor efficiency results in a 2-percent increase in engine 
specific fuel consumption (SFC) and a 2.8-percent loss in engine power at design speed. 
These losses become more pronounced at part-power settings. In addition to good per- 
formance at off -design flow rates it is important that the compressor operate stably over 



the range of flows and speeds required by the engine operating envelope. The usable 
range of the compressor pressure ratio-mass flow characteristic is bounded by the surge 
and choke mass flow rates. Operation at flows less than the surge point flow should be 
avoided because of potentially dangerous vibrations induced by the intermittent flow re- 
versals and power loss. Operation with the compressor choked is generally avoided be- 
cause of the poor compressor efficiency and pressure ratio at the choke point. The 
problem undertaken in this analysis is to determine the centrifugal compressor perform- 
ance characteristics over a range of rotative speeds and flow rates and predict the 
usable range of flow rates at which the compressor can operate. The method of analysis 
uses the loss correlations and equations of reference 1 with surge and choke criteria 
added to predict compressor operating range. 

A FORTRAN program has been developed at NASA Lewis Research Center which 
predicts centrifugal compressor performance through utilization of empirical correla- 
tions which are related to the compressor geometry and velocity diagram characteris- 
tics. A complete knowledge of the compressor overall geometry and working fluid inlet 
total conditions is required for its use. Working fluid state conditions and flow proper- 
ties are calculated using a mean streamline one -dimensional analysis. The program is 
limited to centrifugal compressors with channel diffusers operating up to their choke 
point. A comparison of calculated and experimental performance is given to demon- 
strate the correlation with limited experimental data. The experimental data presented 
for comparison were obtained from a compressor developed by Solar Aircraft Company, 
Division of International Harvester, Inc. Shock losses in the rotor are neglected. 
Clearance losses are considered to be inherent in the impeller losses since good per- 
formance correlation is achieved on compressors operating with reasonable clearances. 


METHOD OF ANALYSIS 

Input values of (V/V ci> )q are used to determine compressor mass flow rate. (Sym- 
bols are defined in appendix A and compressor stations are shown in fig. 1. ) Individual 
losses are calculated using velocity diagram characteristics and empirical correlations 
determined by the input absolute velocity level and compressor geometry. For com- 
pressors having inlet guide vanes it is assumed that the vanes are placed in a constant- 
area annulus having no wall curvatures or slopes. Overall compressor efficiency, total 
pressure ratio, and mass flow rate are tabulated for each operating point inside the pre- 
dicted range for each speed line that is input. A complete discussion of the loss calcula- 
tions is discussed in appendix B. 
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Figure 1. - Meridional cross section of compressor. 


Required Information 

The compressor data needed to use the FORTRAN program are the following: 
inducer inlet tip diameter, inducer inlet hub-tip diameter ratio, inducer inlet blade angle 
at the root -mean-square diameter, inducer blade blockage at the rms diameter, inducer 
inlet hub and tip wall slopes, inducer inlet hub and tip wall curvatures, inlet solid -body 
swirl angle at the rms diameter, inducer -tip to impeller -exit diameter ratio, impeller 
exit blade angle, impeller blade exit height, number of impeller blades at exit, vaneless 
diffuser diameter ratio, vaned diffuser setting angle, total vaned diffuser throat area, 
vaned diffuser area ratio, skin friction coefficient of the wetted surfaces, design rotative 
speed, and the fraction of design rotative speed for which the calculations are to be per- 
formed. 

The working fluid properties and state conditions required are: inlet total pressure, 
inlet total temperature, inlet dynamic viscosity based on total conditions, specific heat 
ratio, and gas constant. 

Two control variables are used to select the options of solid-body prewhirl and im- 
pellers with splitter blades. Inducer inlet hub and tip wall slopes and curvatures are all 
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set equal to zero for compressors with inlet guide vanes. 

Information which corresponds to the design speed line for the example compressor 
is as follows: 


Compressor data 


Design rotative speed, N, rpm 72 000 

Inducer inlet tip diameter, D^, m 0.0813 

Inducer hub-tip diameter ratio, X 0.5313 

j Inducer inlet hub wall curvature, C^, m~^ (positive around axis of rotation) 0 

Inducer inlet tip wall curvature, C 1T , m' 1 (positive around axis of rotation) 0 

Inducer inlet hub wall slope, Xjj, deg from axial 7 

Inducer inlet tip wall slope, x^, deg from axial 0 

Inducer inlet blade angle at the rms diameter, deg from meridional 49.0 

Inducer -tip to impeller -exit diameter ratio, D 1T /D 3 0. 5614 

Number of impeller blades at exit, 32 

Impeller blade exit angle, deg from meridional (positive opposite to direction of rotation) 25 

Impeller exit blade height, b^, m 0.0051 

Vanless diffuser diameter ratio, D^/D^ 1.14 

Vaned diffuser setting angle, a A , deg from meridional 78 

q 2 

Total vaned diffuser throat area, A^,m 0.00071 

Vaned diffuser area ratio, Ag/A^ 2.7 

Friction coefficient, 0.004 


Working fluid properties and state conditions 

Specific heat ratio, y 

Gas constant, R, J/(kg)(K) 

o 

Inlet total pressure, p^, N/m 

Inlet total temperature, T' K 

u 2 

Inlet dynamic viscosity based on total conditions, N-sec/m 

Velocity diagram properties 


Number of values of compressor inlet absolute critical velocity ratios used in calculating performance . . 15 


\ ci 7q 

Inducer inlet solid-body swirl angle at the rms diameter, deg from meridional 0 

Splitter control variable, SPLT (SPLT = 0, no splitter blades; SPLT = 1, splitters) 1 


. . . . 1.4 
. . 287.05 

. . 101 325 
. . 288.15 

1.788X10' 5 


4 







Calculated Losses 


Inlet guide vane loss . - Equations relating working fluid kinetic energy level, bound- 
ary layer, and blade geometry were developed for axial turbine stators in reference 2. 
These equations were used to calculate state conditions downstream of the inlet guide 
vanes; thus the velocity diagram characteristics were found by satisfying one- 
dimensional continuity . 

An inlet guide vane loss coefficient is calculated using the same assumptions about 
blade geometry and boundary layer characteristics as used in reference 3. The inlet 
guide vane loss is then calculated as a fraction of the ideal kinetic energy at the rms 
diameter. 

Inducer incidence loss . - Reference 4 developed equations for optimum inducer in- 
cidence angle of centrifugal compressors with axial inlets. The optimum incidence angle 
is calculated from the inlet velocity diagram characteristics and blade blockage at the 
rms diameter assuming an incompressible working fluid. The enthalpy loss due to inci- 
dence is calculated assuming that the relative velocity component normal to the optimum 
incidence angle is lost. Therefore, the incidence loss is simply 


Ah INC ” 



where 


W L “ W 1MF sm l^opt " ^lMF I 

Blade loading loss . - Boundary layer growth in the impeller is highly dependent on 
the diffusion of the working fluid internal to the impeller itself. Reference 5 proposed 
an equation for calculating the diffusion factor of the impeller based on a uniform veloc- 
ity loading along the blade chord. This equation is used to calculate the impeller dif- 
fusion factor for impellers without splitters. A modified form, with reduced penalty 
due to aerodynamic work input, is used for impellers which have a set of splitter blades. 
With the diffusion factor calculated by these methods the blade loading loss was ex- 
pressed as 


ah BL “O ' 05 DfUj 

Skin friction loss . - In addition to the losses resulting from the aerodynamic load- 
ing of the impeller blades, the impeller incurs losses due to skin friction of the impeller 
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and shroud wetted areas. Reference 5 developed an equation for this loss based on fully 
developed turbulent pipe flow. In the case of impellers with splitter blade rows, the 
empirical constant appearing in the equation is modified to account for the higher mean 
channel relative velocity caused by the addition of splitters. The general equation used 
for skin friction loss is 


_L_ 

D 0 


Ah SF " K DF C f 


D 


HYD 


D, 



u 


2 

3 


where Kg F =5.6 for conventional impellers and Kg F = 7.0 for impellers with tandem 
blades. 

Disk friction loss . - The specific loss due to windage on the compressor back face 
is calculated using the equation 


Ah DF = 0.01356 


P3 


w Re 


0.2 


u 3 D 3 


This is a form of the disk friction power loss of reference 6. 

Recirculation loss . - Losses resulting from work done on the working fluid due to 
backflow into the impeller are expressed as 


Ah RC = 0.02 j/tan <* 3 

This is a modification of the equation proposed in reference 5. 

Vaneless diffuser loss . - The flow angle and Mach number variation with radius in 
the vaneless space are determined by numerical solution of the differential equations for 
vaneless space flow developed in reference 7. The equations are simplified through the 
assumptions of adiabatic flow and constant geometric depth passage. The radial total 
pressure distribution in the vaneless space is calculated using the equation derived in 
reference 5. When the fluid state and flow properties are determined at the vaned dif- 
fuser leading edge radius by the methods described previously, the vaneless diffuser 
loss is calculated from the equation 
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Vaned diffuser loss . - Curves of maximum static pressure recovery coefficient at a 
given area ratio were extrapolated from the test data for square throat diffusers reported 
in reference 8 for various combinations of vaned diffuser throat Mach number, aerody- 
namic blockage, and area ratio. These data were recorded for channel diffusers with 
symmetrical pressure loadings about the channel centerline. The vaned diffuser in a 
centrifugal compressor is loaded with a pressure gradient across the channel. In this 
analysis it was assumed that the difference in the loadings between the test diffusers of 
reference 8 and the compressor diffusers would have no effect on calculated diffuser re- 
covery. A one -seventh power velocity distribution in the boundary layer along the vane- 

less space endwalls is used to calculate the displacement thickness representing vaned 

* £ 

diffuser throat blockage. The value of static pressure recovery coefficient Cp cor- 
responding to the vaned diffuser geometric area ratio, inlet Mach number, and aerody- 
namic blockage is extrapolated from the test data. 

The vaned diffuser exit critical velocity ratio is calculated using one -dimensional 
continuity . Vaned diffuser loss is then calculated using the equation 


Ah VD " C p T 3 




Calculation of Compressor Choking Flow 

Two criteria are used in the prediction of compressor choking flow: (1) inducer 
choke and (2) vaned diffuser choke. Reference 9 presented data on inducer choking inci- 
dence levels as functions of relative Mach number and blade angle at the rms diameter 
for inducer pitch-chord ratios of 0. 5 and thickness -chord ratios of 0.05. These data 
were used to predict compressor choking flow due to inducer choke. 

Vaned diffuser choke is predicted from one -dimensional continuity using the com- 
puted values of weight flow, total temperature, total pressure, aerodynamic blockage, 
and the geometric throat area. The maximum value of the one -dimensional weight flow 
function is calculated from the equation 
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At the first input critical velocity ratio where the left side of the equation equals or ex- 
ceeds the right side the compressor is choked. One would expect the speed lines to 
have infinite slopes at the choke point. The vertical portions of the speed lines are ap- 
proximated through the assumption of zero static pressure recovery in the vaned diffu- 
ser for the cases of vaned diffuser choke. The impeller characteristic should exhibit 
a vertical slope at the impeller choke point but this is not predicted by the correlations 
used. 


Calculation of Compressor Surge Point 

The vaned diffuser is assumed to be the component which governs the location of the 
compressor surge point. Examination of performance characteristics of centrifugal 
compressors covering a range of pressure ratios indicates that compressor flow range 
can be expressed as a function of the vaned diffuser leading edge Mach number. With a 
good impeller -diffuser throat area match the flow range will vary with the diffuser lead- 
ing edge Mach number as shown in figure 2. This curve was deduced from unpublished 
operating range data for centrifugal compressors with design point pressure ratios 
ranging from 1.9 to 10.0. 



Figure 2. - Compressor flow range as a function of vaned diffuser leading edge 
Mach number. 
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Comparison of Estimated and Experimental Performance 


This section compares the results obtained by the analytical procedure described 
herein to the experimental results obtained with a commercially manufactured centrifu- 
gal compressor. 

Figure 3 shows the comparison between calculated and experimental pressure ratios 
as functions of corrected flow for 60, 80, 90, and 100 percent rotative speeds. Choking 
flow is caused by inducer choke at 100 percent speed and by diffuser choke at the other 
speeds. The assumption of zero static pressure recovery at the vaned diffuser choke 
point permits an estimate of the flow rate at which the compressor pressure ratio char- 
acteristic becomes vertical. Experimentally, the impeller characteristic exhibits be- 
havior similar to that of the diffuser at the choked condition. However, the multi- 
dimensional effects that characterize this behavior apparently cannot be approximated 
accurately with strictly one -dimensional correlations. Modifying the correlations for 
evaluation at the hub, mean, and tip sections and averaging the results may help to pre- 
dict the gradual pressure ratio reduction in moving toward choke. The potential refine- 
ments offered by this approach have not yet been considered herein. Maximum error in 
predicting the choking weight flow occurs at the 90 percent speed line and is 25 percent. 
Maximum error in predicting the compressor operating range ( w sur g e / w c hoke^ occurs a ^ 
the 80 percent speed line and is 8. 1 percent. 

Figure 4 shows the comparison between calculated and experimental efficiencies as 



Figure 3. - Comparison of calculated and experimental 
compressor total pressure ratios as functions of equiv- 
alent weight flow and rotative speeds. 



Figure 4. - Comparison of calculated and experimental com- 
pressor total efficiencies as functions of equivalent weight 
flow and rotative speeds. 
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functions of corrected flow. The worst estimate of peak efficiency occurs at the 90 per- 
cent speed line and is in error by 1. 2 points. Only partial data is shown for the effi- 
ciency contour at 60 percent speed. This results because the data is replotted from a 
performance map that had efficiency contours and more data could not be accurately in- 
terpolated. Predicted surge point efficiencies for all speed lines are considerably lower 
than those experimentally measured. This can be attributed to predicted diffuser pres- 
sure recovery being lower than that attainable experimentally. Aerodynamic blockage 
growth may not be as rapid as that predicted by the model in moving toward surge. 

Compressor efficiency increases with increasing flow until the choke point is 
reached. At the choke point vaned diffuser static pressure recovery decreases, result- 
ing in a decrease in overall compressor efficiency. At the impeller choke point one 
would expect the impeller losses to increase substantially, resulting in a decrease in 
overall compressor efficiency. However, this behavior is not predicted through the use 
of the empirical correlations. 


CONCLUDING REMARKS 

This section summarizes the pertinent characteristics of the prediction method, 
discusses the limitations of the method, and notes the effect of the limitations on the re- 
sults. 

For inducers with relative Mach numbers appreciably in excess of unity the pre- 
dicted compressor total pressure ratio and efficiency are greater than that attainable 
experimentally. This is a result of neglecting shock losses in the inducer inlet. 

Predicted static pressure recovery for vaned diffusers with throat Mach numbers 
greater than unity will be too large . Extrapolation of the channel diffuser data used in 
the program is valid only in the subsonic flow regime. 

The prewhirl option available in the program is for solid-body vortex swirl. This 
type of swirl is used primarily because it minimizes inducer inlet tip relative Mach 
number . 

Compressor total pressure ratio and efficiency are not adequately predicted with 
the one -dimensional correlations when the compressor mass flow rate is limited by the 
impeller rather than the vaned diffuser. The predicted total pressure ratio and effi- 
ciency levels are higher than those measured experimentally. 

Predicted surge point efficiencies are considerably lower than those measured on 
the example compressor. Correlation with experimental data from other compressors 
indicates that the predicted surge point efficiencies are sometimes higher and some- 
times lower than the corresponding experimental values. This may be due in part to the 
difficulty encountered in making precise measurements at this flow condition. 
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FORTRAN PROGRAM 


Program Input 


The information required for the program is input using NAMELIST format. A 
sample input sheet, with the data used for generating the 100 percent speed line for the 
example compressor, is shown in figure 5. The definitions of the FORTRAN input varia- 
bles are: 


ADTH 

AL1MF 

AL3 

AR 

B1MFB 

B2 

B2X 

BLOCK 


2 

total vaned diffuser throat area, Ag, m 

inducer inlet absolute flow angle for solid-body swirl at rms diameter, 
a lMF 

vaned diffuser setting angle, a deg from meridional 
vaned diffuser area ratio, Ag/A^ 

inducer rms blade angle, deg from meridional 

impeller exit blade height, bg, m 

impeller blade exit backsweep, /3g^, deg from meridional 
inducer rms blade blockage function 
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CF 

CHIH 

CHIT 

CURVH 

CURVT 

DRAT 

DIT 

GAM 

LA MX 

MUO 

N 

NONDES 

NVOVCR 

POP 

RGAS 

SPLT 

TOP 

VLDRR 

VOVCR 

Z 


skin friction coefficient, Cj 

inducer inlet hub wall slope, Xjjj, deg from axial 

inducer inlet tip wall slope , x^ , deg from axial 

inducer inlet hub curvature, Cjjj, m”* 

inducer inlet tip curvature, C^rp, m - * 

inducer-tip to impeller-exit diameter ratio, D^rp/D^ 

inducer tip diameter , D^, m 

specific heat ratio, y 

inducer hub-tip diameter ratio, X 

n 

total inlet dynamic viscosity, Pq, N-sec/m' 5 
compressor design rotative speed, N DE g, rpm 
percent compressor design rotative speed, N/N DE g 
number of input values of (V/V )q 

n 

total inlet pressure, Pq, N/m 
gas constant, R, J/(kg)(K) 

control variable for splitter calculations (SPLT = 0, conventional 
impeller; SPLT = 1, impeller with splitters) 

inlet total temperature, Tq, K 
vaneless diffuser diameter ratio, D^/D^ 
values of inlet critical velocity ratio, (V/V cr )g 
number of impeller blades at exit, 


Sample Output 

In the normal mode of operation the program performs two iterations over the 
range of inlet absolute critical velocity ratios which are specified. On the first iteration 
the compressor choking flow rate is determined. On the second iteration the printout of 
compressor performance at flow rates outside the usable range is deleted. A listing of 
the input information is tabulated for the compressor configuration studied. Then the 
compressor surge and choke flow rates are printed out. Finally, for each calculated 
point inside the operating range, compressor equivalent weight flow, total pressure 
ratio, total efficiency, and individual efficiency decrements are printed out. Sample 
output for the example compressor is as follows: 
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Definitions of the abbreviated FORTRAN output variables are: 

DETABL efficiency decrement due to impeller blade aerodynamic loading, 
DETADF efficiency decrement due to disk friction, 

DETAIGV efficiency decrement due to inlet guide vane loss, 

DETAINC efficiency decrement due to inducer incidence, 

DETARC efficiency decrement due to recirculation, 

DETASF efficiency decrement due to skin friction, 

DETAVD efficiency decrement due to vaned diffuser losses, Arjyjj 
DETAVLD efficiency decrement due to vaneless diffuser losses, A^yjj) 
ETAT compressor total efficiency, r; 

WEQ equivalent flow rate, w |/iy6, kg/sec 


Error Messages 

If the compressor choking flow rate is not reached with the specified values of inlet 
absolute critical velocity ratios, the operating range cannot be determined with the 
methods used in this analysis. When this happens, the calculated performance output 
is deleted and the message COMPRESSOR CHOKING FLOW HAS NOT BEEN REACHED 
is printed. 

If the first input value of results in calculated weight flow greater than 

the compressor choking weight flow, the message VOVCR ARRAY TOO LARGE is 
printed . 

For input values of ^V/V cr ^ much too small for the compressor of interest a 

velocity diagram solution at the impeller exit cannot be obtained. When this occurs 
the message IRRATIONAL EXIT TRIANGLE is printed. All of these problems can be 
corrected by making the appropriate adjustments to the VOVCR array used. 


Main Program FORTRAN Variables and Corresponding Engineering Symbols 

A1 a^ 

A2 a 3 

ACOUSR a /a' 
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ADTH 

AL1 

AL1MF 

AL2 

AL3 

ALPH1 

ALPHA 

AMT 

AMU 

A MULT 

ANU 

AR 

B 

BO 

B1 

B1H 

B1MF 

B1MFB 

B2 

B2X 

BARR 

BLOCK 

BMULT 

BOPT 

BT 

CF 


“IT 

“IMF 


“4 

temporary storage 
a 



Ht 

Hu 

^1MF 

^lMFb 

b 3 

^3b 

array of vaned diffuser throat blockages 



^lMFopt 

vaned diffuser inlet blockage 
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CHIH 

X 1H 

CHIMF 

X 1MF 

CHIT 

X 1T 

CINC 

inducer choking incidence 

CINC1 

interpolation variable 

CINC 2 

interpolation variable 

CINC40 

choking incidence array for 40° rms inducer blade angles 

CINC 50 

choking incidence array for 50° rms inducer blade angles 

CINC60 

choking incidence array for 60° rms inducer blade angles 

CONST 1 

k bl 

CONST 2 

CO 

COSA 

cos a ^rp 

CP 

s 

CPSTAR 

p* * 

^p 

CURVH 

C 1H 

CURVMF 

C 1MF 

CURVT 

C 1T 

D1H 

D 1H 

D1T 

d it 

D2 

°3 

DEBL 

A? ?BL 

DEDF 

«P 

< 

DEEXIT 

A?7 exit 

DEIGV 

At ?igv 

DEINC 

At? inc 

DELTA 

6 

DELTAB 

AB 

DELTAR 

Ar 

DELTAS 

A6* 
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DERC 

At? rc 

DESF 

A7? sf 

DEVD 

At? vd 

DEVLD 

At? vld 

DF 

D f 

DHACT 

Ah act 

DHAERO 

A ^aero 

DHBL 

Ah BL 

DHDF 

Ah DF 

DHDIF 

Ah VD 

DHEST 

Ah est 

DHIGV 

Ah IGV 

DHINC 

Ah INC 

DHRC 

Ah RC 

DHSF 

Ah SF 

DHVLD 

Ah VLD 

DHYD 

D HYD / ^ D 3 

DIFLEM 

array 

DMF 

°1MF 

DRAT 

D 1t/°3 

EPS 

e 

ETAD 

^AD 

ETAR 

^R 

EWF 

w ^0^6 

FI 

C* * at vaned 

F2 

C** at vaned 
p 

F3 

C* * at vaned 

F4 

C* * at vaned 

F5 

C* * at vaned 
P 


diffuser area ratio of 1 . 2 
diffuser area ratio of 2.0 
diffuser area ratio of 3.0 
diffuser area ratio of 4.0 
diffuser area ratio of 5.0 
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FINT 


C dn 


FLFUNC 

FLRNG 

G1 

G2 

GAM 

GINT 

HO 

HI 

HIS 

INC 

J 

JAN 

K 

L 

LA MX 

LOD 

MIN 

MUO 

MWRMS 

N 

NO 

NONDES 

NVOVCR 

NWRITE 



r lMF 

1H 




(r+l)/(r-l) 


w surge //w choke 

y + 1 

r - l 

r 



C dn 


r lMF " r lH 
r lT " r lMF 
Ah is 

%MF " ^lMFb 


index variable 
iteration counter 
choke indication variable 
index variable 


X 

L/D 3 

index variable 
^0 

array of inducer rms relative Mach numbers 
N DES 

index variable 


N/>N DES 

number of values of (v/V cr ^ 
output control variable 


0 
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POP 

P 0 


PI 

P 1MF 


PIP 

P 1MF 


P2 

p 3 


P2P 

p 3 


P3P 

p 4 


P4 

p 6 


P4PG 

iteration variable 

P4PL 

iteration variable 

PEXIT 

p 6 


PHI 

<P 


PIMF 

P 2MF 


PIPMF 

P 2MF 


POPP1 

( p / p ')lMF 

PP 

storage variable 

PPEXIT 

p 6 


PPOP1 

(P'/pliMp 

PR 

p 6 /p 0 


PREC1 

array of 

C!** at area ratio 1. 2 

PREC2 

array of 

C** at area ratio 2.0 

PREC3 

array of 

C* * at area ratio 3.0 

PREC4 

array of 

C** at area ratio 4.0 

PREC5 

array of 

C** at area ratio 5.0 

PTH 

p 5 


PTHP 

p 5 


Q 

Q 


QAERO 

a 

^aero 


R 

r 


ROP 

p 0 


R1 

p l 
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RIP 

P 1 

R2 

r 3 

R2G 

P3,est 

RANGE 

w surge //w choke 

RARRAY 

temporary storage 

RE 

Re 

RGAS 

R 

RH02 

P 3 

RHOH 

PlH 

RHOMF 

PlMF 

RHOR 

p/p' 

RHOT 

pit 

RMRMS 

w imf/ & i 

RTHP 

Ph 

RVOEST 

KP' V cr) 2 

RVORP1 

(pV/p’Vc^ 

RVORVC 

(pV/p-V,^ Pi/P 1MF 

RVRVC4 

(pv/p-v cr ) 6 

RVRVTH 

(P v /P’ v cr) 5 

S 

9 

integral of M 

SINA 

sin 

SPLT 

splitter control variable 

SW 

w 

TO 

T 0 

TOP 

T o 

TOPP 

T! 1 

x o 

T1 

T 1 

T1PP 

rp f f 
1 
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T2 

T2P 

T2PEST 

T2PP 

TALPH 

THETA 

TOTO 

TOT1 

TPL 

TX 

U1H 

U1MF 

U1T 

U2 

V1MF 

V1T 

V2 

VARAL 

VARAL1 

VARM 

VARM1 

VCR 

VCRTH 

VLDRR 

VM1 

VM1H 

VM1MF 

VM2 

VM1MF 


T 3 

T 3est 

Tt f 

3 

tan a 

e 

T 1 /T 0 

V T 1 

»'/»3 

"R ih aero/ C p T o * 1 


U 


u 


u 


1H 

IMF 

IT 


u. 


IMF 

r 

IT 


A tan a 

temporary storage 
AM 2 

temporary storage 

(Mo 

( v 4 

04/03 

V mlT 

V mlH 

V mlMF 

V m3 

V mlMF 
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VOVCR 

( V/V cr) Q 

VOVCR4 

( V/V cr ) 6 

VOVCR1 

^cr)j 

VSL 

V SL 

VU1H 

V ulH 

VulMF 

V ulMF 

VU1T 

V ulT 

VU2 

V u3 

W1H 

W 1H 

W1MF 

W 1MF 

W1MFEF 

W lMFeff 

WIT 

w iT 

W2 

W 3 

W20W1T 

Wg/WiT 

WCHK 

w choke 

WCR 

W cr 

WFUNC 

w V F 3/ t 5 A 

W1MF 

W 1MF 

WL 

W L 

WOU2 

(W/Ug ) 2 


\ J/ av 

WSURGE 

W 

surge 

WU1H 

W ulH 

WU1MF 

W ulMF 

WU1T 

W ulT 

WU2 

W u3 

XI 

oi - a 4 
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xc 

K vu 

XK 

K vm 

XM 

M 

XM1 

temporary storage 

XM2 

m 3 

XMACH 

M 5 

XMARR 

array of vaneless diffuser Mach numbers 

Z 

Z 3 

ZETA 



Program Listing 


DIMENSION VOVCR (15) 

OIMENS ICN P3P4 20) ,XMARRI20> 

DIMENSION E ( 20 ) » S 4 20) 

OIMENS ION AMT44),8ARRl6)«PRECi44,6), PREC2 14*6) ,PREC3(4,6) 

DIMENSION PREC444.6) ,PREC5<4,6) 

D I MENS ION MtaRMS48).ClNC6048> .CINC5018) ,CINC4 0< 8) 

DIMENSION FLRNG4 7 )«DIELEM4 7) 

INTEGER SPLT 

REAL LAMX , LCD. MUO.NUNDES.MWRMS , INC .N 
OATA4AMT4 [)«1 3 1«4)/»2i*^9»6i*8/ 

DATA! BARR 4 l ). l=1.6)/.02. ,04,.06,.08,. 10, . 12/ 

DATA! 4PREC14 I . J), 1=1 • 4 ) « J» 1,6 ) /.234, .244, . 257* .269 • .215, .224, .233, 

1.243. . 207. .215. .223. .232. . 193 .. 199 . .206 . .2 12 .. 183. .190 . . 196 . .202 . . 

2169. . 1 76. .182. .188/ 

DATA! I PR EC 2 I I . J ) . I=1,4),J=1,6) /.644, . 67 J , .696, .722 , .620, .638, .656, 
1.674, .590, .606, .623. .639, . 562 ,. 976. . 590, .605 ,. 538 , .551 , . 564 , . 578 , . 
25 1C. .9 24, .538, .552/ 

DATA! I PREC3 ( I , J ) . 1*1 ,4) * J*1 ,6) /. 782, . 789, • 796, .802 ,. 750, .756, .762, 

1.766. . 7C8 « .716, .724**7 32 , . 67 2,. 680, .687, .695 , . 652 , .648 , .654 , .660 , . 

2604.. 6 12. .619. .626/ 

DATA! ( PREC4I I. J ), 1= 1 . 4 ) , J= 1 ,6 ) / . 842. . 838, . 833. .828 , . 8. .8 , . 8 , .8 , .75 

12. . 7 56. .7 60. .763. .7 10. .7 13. .716. .7 19 . . 675 . .678 . .680. .683. .630. .635 

2. . 640. .646/ 

DATA! 4PREC54 I.J), 1=1,4). J= 1,6)/. 878. .865,. 852, .838 ,. 832 , .825 , .818, 
I • 6 1 2, . 78, . 78, .7 8, .78, . 736, • 735 , .735, . 7 34 , .692, .694 ,.695, .696, .644, 
2.647. .650. .652/ 

D A T A I 0 I NC6 0< 1). I — 1,8)/— 10. 5, -5. »-2 ., 1.5, 2. 5, 4. ,5. ,5./ 

0 A T A 4 C I NC 504 I) , I — l,6)/~15., w 7.5, v 3*,l«5.4.,6.5,9«,ll./ 

OAT A( C INC40 4 I ) . 1= 1. 8 ) /-20 . .- 1 2 . ,-5 . , 1 . 5, 5. , 8 . , 11 .5 , 14. 5/ 

OAT A4 MtaHMS 4 I). 1 = 1, 8)/. 5, .6, .7, .8, .9,1. ,1.1, 1.2/ 

OAT A 4 D IELEM4 I). 1=1, 7)/ .4, .5,. 6. .7, .8, .9,1./ 

OAT A! ELRNG ( I), 1 = 1, 7)/. 48, .56, .64, . 72 . .79. .84 ,. 89/ 

NAMELIST /INPUT/ GAM , P OP . TOP . N ,U1T ,MUO »CE . WO VCR, NVOVCR , RGAS * 

10k AT ,1 AMX.82X.Z .VLDRR. 82, B1MEB. AR. BLOCK. AL3. AOTH,NONDES, SPLT, 

2 AL IMF. CHIB. CHIT. CURVH ,CUR VT 
1 RE AC 4 9 , INPUT ) 
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MR l TE I 6 * 309 * 

MRITEI6.315) P OF. TOP .MUO.GAM.RGAS 
mR ITE (6*3 10) 

mR ITF<6.316> l .CF.SPLT.N.AL IMF, BLOCK 
MHITE<6.311 I 

MR ITE(6,317) B1MFB.LAMX.0RAT 
MR ITFI6.312 ) 

MRITFI6.318) B2X.B2. VL0RR.AL3 
MR I TF 1 3 13 ) 
mRITE<£.319> AOTH.AR 
MR ITE 4 6.3 14 > 

MR ITE <6 .320) ( VOVCR ( 11.1=1. NVO VCR) 

PEROS=NCNDES*100. 

G1=GAM*1. 

G2=GAM— 1 • 

CP =GAM*RGAS/G2 
H2X=.01745*B2X 
AL 3=AL3*.01745 
AL 1MF=.01745*AL IMF 
CHI H= .0 1745*CH I H 
CHIT=.C1745*CH£T 

FLFUNC=S0RTIGAM/RGAS*(2. /G1 )**(G1/G2I > 

ROP=POP/RGAS/TOP 
VC R= SORT ( 2 •*GAM/G1*RGAS*T OP) 

UIT=3. 14159*N*N0N0ES*01T/60. 

NMR 1TF=1 
00 18 L = 1 • 2 
M I N= 1 

2 OC 17 J = M IN .NVO VCR 
K=C 

BIMF8=.01745*B1MFB 
VO IT=0 • 

VU1MF=0. 

VU1H=0. 

VM l H= VO VCR I J )*VCR 
U2=G1T /DRAT 
02=D1T/CRAT 
O 18=0 1 T*L AMX 
018=0 1 T*L AMX 

OMF=SORT I D 1T**2*I 1 .«-LAMX**2 » /2 . * 

01MF=U1T*GMF/0 IT 

CH IMF= (CHI T-CH I H ) * I DMF— 0 1H ) / ( 0 IT— 0 1H 1 +CH I H 
COR VMF = ( COR VT— CUR VH ) * I DMF-018 1 / I D 1 T-0 1H 1 +CUR VH 
H0= I CMF— D 1 F ) / 2 • 

Hl=IClT-DMF>/2. 

F IKT=HO/2.*(CURV8+CORVMF > 

AMMT=EXP(FINT) 

vmimf=arolt*vmih 

GIM = I H0+H1 1 /B.*( I 2.— HI /HO 1 *CUR VH+( HQ+H 1 1**2 /HO/Hl*CURVMF+( 2 «-H0 

1/81 >*CORVT> 

BRLLT=FXP(GINT> 

VM1T=BH0LT*VM1H 

VMlhN=VMlH*COS(CHlH) 

VMIMFN = VM1MF*C0S< CHlMF) 

VM1TN=VM1T*C0S(C8IT) 

R8O8=K0P*< 1 .—VM1H** 2/2. /CP/TOP >** I 1./G2) 

K8CMF=R0P*( l.-VMlMF**2/2./CP/T0P)**( 1./G2I 
R8GT=R0P* I l.-VMlT**2/2./CP/TCP 1**1 1./G21 
FCM=KFG8*VMlHN*ClH/2. 

FCK2=R8CMF*VMlMFN*DMF/2. 


27 



FCN 3=RFCT*VRlTN*01T/2. 

Sta= ( H0 + H1 ) /6.*4 (2.-Hl/H0>*FCNl+<H0+Hl >**2/H0 /H 1*FCN2+ 4 2 . -HO/H 1 >* 

1FCN3)*6. 28318 
DF IGV=0 • 

V 1 RF = SORT ( VH IMF **2+ VU 1 MF**2 > 

Tl=T0P-VlMF**2/2./CP 

P1P=P0P 

POPP 1=4 l.-VlMF**2/2./CP/TOP»**4GAM/G2> 

P 1=P1P*P0PP1 

R1F=PIP/RGAS/T0P 
R I * R1P*4P1/PIP>**< i./GAM) 

RE=U2*Q2/MU0*R0P 
A1 STAG=AL IMF/2. 

F S=0 .0076 / 4 COS 4 AL lMFl-.025>*4 1 .♦C0S4 ALSTAGJ/0.7) 

V0VCRl=V0VCR4 J) 

IFIALIRF.LT. 0.001) GO TO 36 
9 VOWCR l = VO VCR 1 *.00 1 
AkE=( VCVCR l*VCRI**2/2. 

POPP 1= 4 l.-AKE /CP/TOP 1**4 GAM/G2) 

AKFIC=AKE/4 l.-ES) 

P10PP0=4 l.-AKEID/CP/TQP)**4GAM/G2l 
P1P=POP*PIOPPO/POPP1 
R 1 F =P IP /R GAS/TOP 
P1=P1P*P0PP1 

R1=R1P*4P0PP1>**4 l./GAM) 

01=3. 14159*01T**2*4 l .-LAMX**2 I *VGVCR l*VCR*C0S4 ALIMF1M. 
t.TCHK = Gl*R 1 

IFIFTCHK.LT. SM> GO TO 9 
tf U 1MF= VGVCR 1*VCR*SIN4ALIMF) 

VM1MF=VGVCR1*VCR*C0S4 AL1NF) 

V1RF=S0RT4 VUIMF**2+VM1MF**2» 

T 1=T0P— VlMF**2/2. /CP 

XK=VU1MF**2*2.*VR1RF**2 

XC=VU1RF/DMF*2. 

Vl)lT = XC*01T/2. 

VUlH=XC*01F/2. 

VR1T=SCRT4 XK-2.*VUIT**2> 

VR1H=S0RT4 XK-2.AVU1H**2> 

0F I GV=ES*AKE ID 
36 VIT = SURT4 VR 1T»*2+VU1T**2 > 

AL 1=ARSIN< VL1T/V1T) 

W01T=U1T-VU1T 

H1 = ATAMFLI1T/VM1T> 

W 1T=SURT4 VM1T**2+HUIT**2I 
RlilH=(JlH—VUlH 
B 1 H = A T AN4 FUIH/VMIH) 
w 1H=SQRT ( VHlH**2*FUlh**2> 

A1=SCRT4GAM*RGAS*T1» 

RU1RF=U1RF-VU1MF 
W1RF=SURT< VM1MF**2+FUIRF**2) 

RRRMS=tolHF/Al 
B1RF=ATAN4 RUIRF/VM1MF) 

INC =4 B1RF-81MFB >* 57. 29 57 7 

£PS=AT AN4 4 1 .-8LQCK)*TAN4 B1MF ) / 4 1.* BLOCK* 4 TAN 4B1MF) »**2) ) 

BGPT=B1RF-EPS 

T1FP=T 1«U lRF**2/2./CP 

WCR = S0RT 4 2 .*GAM/G l*RGAS*T 1PP) 

T4~. T0= I .— G2/Gl*(HlMF/hCRI**2 
WlHFEF = tal)4F*C0S( B0PT-B1MFB ) 
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T0T1=I.-G2/G1*( 8lHFEF/8CRl**2 
T 1 = 1 1PP*T0 T 1 

8L=R1MF*SIM A8SIB0PT-B1MF1 ) 

0HINC=8L**2/2. 

P IFMF=PlP*EXPl-OH INC/T I/RGAS1 
DEL TA=PCP/ 101325.35 
TFETA=T0P/2B8. 15 
E8F=S8*S0RT1 THETA) /DELTA 
81HFB=B1HFB*57.29577 

I F i B1HF8.GE .40. • AND.B1HFB.LE ■ 50. > GO TO 71 

1 = 2 

81 I F ! RHRHS .GE .H8RHS11— 1) • AND .RHRHS. LE. M8R MS! III GO TO 80 
1 = 1 + 1 

GC TO 81 

80 C INC1= ( RHRHS— M8RH SI 1-1 ) 1/ ( H8RHS1 I 1-HWRHS1 1-1 II MCINC60I I 1-CINC60 1 I 
1-1 > 1+CINC601 1-11 

C INC2=!RHRHS— H8RHS1 1-1) )✓ ( H8RHS1 I 1-M8RHS1 1-1 ) ) ♦( C INC501 I >-C INC 501 I 
1-11 )+CINC50( I-ll 

C I NC=C INC 1 + 1 60.-B1HFB) / 1 0.* ( C 1NC2-C INC 1 1 
IF 1 INC.LE.CINCI K = 1 
GO TO 74 

71 1=2 

72 IFIRHRHS.GE.H8RHS1 I-ll .AND. RHRHS. LE.HNRMS1 II 1 GO TO 73 
1=1 + 1 

GC TC 72 

73 C INC 1= ( RHRHS— H8RH SI 1-1 1 1/ (H8RMS1 1 1-H8RHS! £-1 1 l*!CINC50i 1 1-CINC50( l 
1— 1 1 1+C INC501 I-ll 

C INC2=!RHRHS-HHRHS1 1-1 1 1/ 1H8RHSI 1 l-HMRHS! 1-1 1 1*1CINC401 1 1-CINC401 I 
1-11 1+CINC401 I-ll 

C I AC=C INC 1 + 1 50.-B1HF8) /l 0 .*1 C I NC2-C INCH 
IF I INC.LE.CINCI K=1 

74 CCA T I NUE 

IFIK.EO.) • ANC. J «E 0. 1 1 GO TO 202 
IF1K.EQ.1.AND.L.E0.1! GO TO 18 
T2PP=T 1PP+IU2**2-U1MF**2 1/2./CP 
PH I=VM1HE/U2 

EPSLIM=1./EXP(8 • 16*C0SI B2X 1/2 1 
VSL=SORT 1C0S182X1 >*U2/Z**.7 

I F I ( CMF/02 ) .GT.EPSLIH1 VSL=U2*I SORT 1C0SIB2X1 1/2**. 71*! l.-l C DHF/D2- 
1EPSLIM1/I 1.-EPSLIM1 )**3l+U2*l < OMF /02-EPSLI Ml /I 1.— EPS LIM1 1**3 

QHEST=U2**2 

T2PEST=IDHEST/CP/T0P+1 .)*T0P 
R2G=RI*1T2PEST/T0P1**1 1./G21 
4 RHC2=R2G 

VH2=S8 / (3. 141 59*R R02*02*B2 ) 

Vl)2=U2— VH2*TAN1B2X1— VSL 
MIN=J 

IF! VH2.LT.0..AND.L.E0.2) GC TO 200 

IF ( VHP «LT . 0. 1 GO TO 17 

IF! VU2.LT. 0..AN0.L.EQ. 2) GO TO 201 

IF I VU2 .IT . 0. 1 GO TO 17 

8L2=t>2-VU2 

82= SORT !8U2**2+VH2**21 

I F ( ( VSL*COS! B2X 1/821 .GT .1.. AN 0. L.EO. 2) GO TO 202 

IF! (VSL*C0S!B2X 1/821 .GT. 1. 1 GO TO 17 

T2=T2PP-82**2/2./CP 

A2=S0RT ! GAH*RGAS*T21 

82C81T=82/81T 

8002=1 PH I **2+1 0HF/02 >**2 + 82081 T**2*( PH I * *2+0 PAT* *2 1 1 /2. 
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AL2^ATAM VU2/VM21 
V2-SORT4 V02**2*VM2**21 
T2P=T2*V2**2/2./CP 

OF A ERO«CP*TOP* 4 T2P/ T OP— 1 • ) 

OAFRC— CHAFRO/U2**2 

DF = * IT /U 2*4 2/ 3. 14159*U.-C1T/D21*2.*01T/D2 1 

CONST 1*0* 1 5 

IF(SPLT.FC.l) CONST 1=0*6 
0F=l.-y20*lT*CQNSTl*UAERG/CF 
DFBL<=0 .05*DF**2*U2**2 
LO0«4 l.-DMF/. 30481 /C0S4B2X1 22. 

DHCF = *0 1356*RH02*U2**3*D2**2/ SN/RE** .2 
DHYD— Z 23 • 14 159/CO SI 82 X) ♦02/ B2 

0FYD*1*/DHYD*01T/D2242*/41 .-LAMX l*2**Z/3* 14159/1 1«*LAMX1 
1*SCRT ( 1**41 .*LAMX**2 1/2.* 4 TAN 4 Bll 1**21 1 
OFRC^O *02* SOR T I TAM AL2 1 1*DF**2 *U2**2 
CONST 2 = 5*6 

IF4SPLT.EC.il CCN ST2*7 .0 

0H$F=CGNST2*CF*LCD/DHYD*NQU2*U2**2 

OHACT=CFAERC*DHDF*CHRC 

H I S*OHAERO— DHBL-DHSF-DHl GV— OH INC 

FTAfi*HIS/OHAERO 

TX = ET AR*CHAERO/CP /TOP* 1. 

51 P2P=TX**4GAH/G21*PIPHF 

P2*P2P*<T2P/T2 1** 4— GAM /G2 1 
R2G=P2/RGAS/T2 

IF4ABS4 (RH02-R2G1/RH021.GT. 0.0011 GO TO 4 
XM2= V 2 / A 2 
R2=02/2 • 

AMU^l .4579E— 6*T2**1. 5/4T2*llC.41 

ANL*AM0/R2G 

BO-1. 

XM*XH2 
ALPFA=AL2 
R— 1.0 

F 4 ll = XM2**3/< l.*G2/2.*XH2**21**4Gi/2./G2l 

P3P I 1 1 = P2P 

XPARRi 1 1-XR2 

OEL TAR = < VLORR-1 . 1 / 10* 

ZETA=CF*R2/B2 

AS-O* 

OC 1000 NQ*2* 10 
XM1 =XH 
ALPH1-ALPHA 

OS-CEL TAR*R2/C0S I ALPHA 1 
AS=AS+OS 

DEL T A S*. 03 7*A S** I -.2 1* < V2/ ANU 1**4 -.21 *DS 

B-B0-2.*0ELTAS/B2 

DEL TAB*80- B 

VARN--2 .* 4 l.*G2/2.*XH**21 /CXH**2-i./COS4 ALPHA1**2) *4 4 GAH*XH**2-T AN 
It ALPHA 1**2 1 ♦/ E T A/ BO/COS ( ALPHA 1*1./ BO*DELTAB/OELTAR-i. /COS 4 ALPHA1 ** 
22/Rl*XM**2*DELTAR 

VAR AL-l./COSI ALPH A 1**2 / 4 XH**2— 41* /COS 4 ALPHA 1 1**2 1* 4 ( 1* *G2 *XM**2 

11*ZFTA/BO/COS4 ALPHA 1*1./ B0*0ELTAB/0ELTAR-XM**2/R1*TAN 4 ALPHA 1*0ELTA 
2R 

VARM1-VARM 
VARAL 1 = V AR AL 
B 1 = B 
B— BO 
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AS* AS— DS 

XM=XM 1**2+ WARM 
XM* SORT 1 XM ) 

TALPH=TAN4 ALPH1 )+ VARAL 1 
AL PH A* A T AMTALPH* 

OS*DELTAR*R2/COS4 ALPHA) 

R=R+OELTAR 

AS=AS+OS 

OFL T AS* «037*AS** 4— .2)*(V2/ANU) **( — »2 ) *OS 
B=80— 2 •*0ELTAS/B2 
DEL TAB = B0-B 

VARM=-2 •* 4 1«+G2/2.*XM**2 ) / 4 XM**2— 1. /COS4 ALPHA) **2) *4 4 GAM*XM**2— T AN A 
14 ALPHA 1**2 )*7ETA/8O/COS4ALPHA)+l./B0*DELTAB/DE LTAR— 1 ./COS I ALPHA)** 
22/R)*XM**2*DELTAR 

VARAL* 1 . /CQS4 ALPHA >**2/ 4 XM**2- ( 1. /COS ( ALPHA) )**2)*( 4 1. +G2 *XM**2 

1**7 FT A /BO /COS ( ALPHA ) + 1 • /BO*DELTAB/DE LT AR-XM**2/R > *TAN ( ALPHA ) *DELT A 
2R 

WARM* 4 WARM 1+V ARM) /2. 

7 WARAL=< WARAL1+VARAL ) /2. 

8=(B+R1 )/2 • 

XM*XMl**2+VARM 
XM*SCRT 4 XM ) 

TALPH*TAN4ALPH1 )♦ VARAL 
ALPHA=ATANCTALPH) 

BO=P 

ACC U SR* SORT 4 1 •/ 4 1 .+G2/2.*XM**2) ) 

RHOR* 1 ./( l.+G2/2.*XM**2)**4 1./G2) 

F4KO)=XM**3*ACGUSR*RHOR*R 

IF4N0.F0.2 > S( NO) =4 F (NO )+F ( NC- 1 ) **0.5*0ELTAR 

IF4NG.EC.2) GO TO 8 

CALL FMTGRL 4 NO. DELTAR.F.S) 

8 TPL=l./4 l.+GAM*CF*R2*S(N0)/C0S(AL2)/B2/XM2*( l.+G2/2.*XM2**2)**( Gl/ 
12./G2) ) 

PP=TPL*P2P 
P3P4NO ) = PP 
XMARR 4 NO ) = XM 
XMACH=XMARR 4 NO ) 

1000 P THP=P 3P 4 NO ) 

1F4L.E0.1) GO TO 30 
M*2 

21 IF4 XMACF.GF.D I FLE M4 M— 1 ) • AND. X MACH. LE.D I FLEM4 M) ) GO TO 22 
M=M ♦ 1 

GC TO 21 

22 RANGE*4XMACH-DIFLEM4 M- 1) ) / 4 01 FLEM4 M)-DIFLEM(M— 1))* 4FLRNG4 M )— FLRNG4 
1M— 1 ) ) +FLRNG 4 M— 1 ) 

WSURGF*RANGE*MCFK 
IF 4 EHF.LT .W SURGE) GO TO 203 
30 CCM INUF 

PTH=PTHP/4 l. + G2/2.*XMALH**2 )**4GAM/G2) 

OHWLD=CP*T 2P*4 4 PTH/PTHP)** 4G2/GAM1-4 PTH/P2P) **4G2/GAM) ) 

BT= 1 •— B 

XMACH=XMACH*C0S4 ABS4 ALPHA-AL3) ) 

PTHP*P TH* 4 1.+G2/2 .*XMACH**2)**4GAM/G2) 

X1=4ALPHA-AL3)*57 .29577 
MFUNC*SM*S0RT4 T2P ) /AUTH/PTHP/B 
IF4HF0NC.GE.FLFUNC) CPSTAR=0. 

IF(toFOAC.GE.FLFUNC) K*1 

IF 4 HFUAC • f E. FLFUNC) GO TO 1* 

IF 4 4 AR-1.2) .GT .0. .AN0.4 AR-2. ) .LT.O. ) GO TO 10 
IF 4 4 AR-2 . ) .GT .0. .AND. 4 AR-3. ) . LT.O. ) GO TO 11 


31 



I F < (AR-3.).GT.0..AND.(AR-4.).LT.0.1 GO TO 12 
If ( ( AR— 4. 1 .GT » 0. »ANO • I AR— 5 • 1 • LT • G • 1 GO TO 13 

10 CALL LIMNT(XMACH,bT,AMT.bARR,PRfcC1.4,6,Fll 
CALL L IMNT ( XMACH » BT • AMT »BARR»PREC2»4#6»F21 
CPSTAR=(AR-1.21/.8*(F2-F11+FI 

GC TO 14 

11 CALL L IMNTIXMACH.bT. AMT, BARR. PREC2. 4.6. F21 
CALL L IN I NT ( XMACH «bT » AMT » BARR * PREC3 »4» 6*F31 
CPSTAR=< AR-2. 1*(F3-F21*F2 

GO TO 14 

12 CALL L IN (NT I XMACH* BT »AMT»6ARR»PRfcC3*4»6»F31 
CA1 L L IN INT ( XMACH *BT tAMT»BARR» PREC4.4# 6 »F41 
CPSTAR=(AR-3. 1*<F4-F31*F3 

GC TO 14 

13 CALL L IN I NT ( XMACH «BT «ANT «6ARR« PREC4.4. 6«F41 
CALL L IMNT(XMACH,BT,AMT.BARR» PREC5.4.6.F51 
CPSTAR=(AR-4.1*(F5-F41*F4 

14 CCNT INUE 

IFlK.EO.l.ANO.J.EC.ll GO TO 202 
IF< K.EO.l.ANC.L.EQ.l 1 GO TO 18 

IFIK. NE.1.ANC.J.E0.NV0VCR1 GO TO 204 
PFX IT=CPSTAR*( PTHP-PTH1+PTH 
RTFiP=PTHP/RGAS/T2P 

VCR TH= SORT ( 2.*GAM/G1*RGAS*T2P 1 
RVRVTH=SM/RTHP/VCRTH/ADTH/B 
VO VCR4= *020 

100 P4PL*PEXIT /< 1 . — G2 /Gl* VO VCR 4**2 >♦*( GAN/G21 

RVRVC4=4 l.-G2/Gl*V0VCR4**21**( 1. /G2 ) *V0 VCR4 
P4PG=R VRVTh*PTHP/ AR/R VRVC4 
VCVCR4=VGVCR4+.00l 
IFIP4PL.LT. P4PG1 GO TO 100 
P PE X IT = P4PG 

IS OhC IF=T2P*CP*( 1 PEXIT/PPEXIT1**(G2/GAM1 — ( PEXI T/PTHP >**( G2/GAM) 1 

PR=PPFXIT /POP 

ET AC=DP AERO-QHSF— CHBL— DFtVLC— OHOI F— DH1GV— OHINC 

ETAO=ETAD/OHACT 

OE I G V=OF IG V 70HACT 

OEBL=DFBL/OHACT 

OESF=OhSF/CHACT 

l)ECF=OFCF / OH ACT 

GERC-=OHRC/OHACT 

l)E VLO=OHVLO/OHACT 

DE VC=OHC 1 F / OH ACT 

OE INC=0HINC/0HACT 

IFIL. FC. 21 GO TO IV 
17 CCNT1N0F 

IFIL.EO.l) GO TO 18 
19 IFINNRITE.NE.il GG TO 35 
MR I T F ( 6 • 300 1 
MRITE46.3011 PEROS 
mR ITE (6.302 1 

MR I TE ( 6 . 303 1 mSURGE.NCHK 
NNR ITF=NmR ITE + 1 
35 M I N=M I N + 1 

MR I TF I 6 .304 1 

mR I TF( 6. 305 1 EHF.PK.ETAO 
MR ITE (A. 3061 

MR ITE (6. 30 71 DE IG V. DE I NC. DfcBL . OESF ,OEDF . DERC .DEVLD .OEVD 
IFlK.NE.l .A NO. MIN.LE.NV0VCR1 GO TO 2 
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IF < K.FC.l .ANC.L .EC.2 ) GO TO 1 
2 00 WRITE16.400) 

MIN=MIN*1 

IF I N IN.GT.NVOVCR) GO TO 1 

GC TO 2 

?Ol HR I TE I 6 .400 ) 

MIN=MIN*1 

IFIMIN.GT.NVOVCR) GO TO 1 

GC TO 2 

202 wRlTE!l.402) 

GO TO 1 

203 WRITEI6.403) VOVCR(J) 

HIK-HIN41 

IFIMIN.GT.NVOVCR) GO TO 1 
GO TO 2 

2 04 HR I TF ( 6*308 > 

GC TO 1 

18 WCHK=SW*SORTI THETA)/DELTA 
GO TO 1 

300 FORMAT ( 59X.13H PERCENT NOES) 

301 FORMAT I63X.F5.1 > 

302 FORMAT i 45X * 16H SURGE FLOW RATE * 10X , 1 5HCH0KE FLOW RATE) 

303 FORMAT I51X.F5.3.20X.F5.3) 

304 FORMAT 1 44X * 4H WEO * 10X. 14HPRESSURE RATIO* 10X, 4HETAT ) 

305 FORMAT I44X»F5.3,13X*F5.3»15X»F4»3) 

306 FORMAT! 1 IX* 8H DETA IGV * 8X • 7H0ETAINC .8X.6H0ETABL *8X*6H0ETASF * 8X.6H0E 
1TACF* 8X «6F0ETARC* 8X* 7H0ETA VLO • 8X« 6HDETAVD) 

307 FORMAT! 1 3X , F6. 5,9X ,F6. 5* 8X ,F6 . 5, 8X *F6. 5, 8X *F6. 5. 8X ,F6. 5. 9X * F6. 5 , 
18X.F6.5.9X.F6.5) 

308 FORMAT ( 45H COMPRESSOR CHOKING FLOW HAS NOT BEEN REACHED) 

309 FCRMATI21H1 INLET TOTAL PRE SSURE , 5X *23H INLE T TOTAL TEMPERATURE , 5X ,2 
19MNLFT TOTAL OYNAMIC VI SCOSI TY.SX .19HSPECIF IC HEAT RAT I 0* 5X, 12HGA 
2S CONSTANT) 

310 FORMAT ! 17H NUMBER OF BLADES , 5X .20HFR ICT I ON CCEFF IC IENT ,5X . 4HSPLT, 5 
IX « 12HDES I GN SPEED*5X.23HINDUCER RMS SWIRL ANGLE , 5X ,27H INDUCER RMS 
2BLCCKAGE FACTOR) 

311 FORMAT ! 1 1 X * 24H INDUCER RMS ELADE ANGLE ,5X . 30HI NDUCER HUB-TIP DIAME 
ITFR RAT 10, 5X.40HINDUCER TIP-IMPELLER EXIT DIAMETER RATIO) 

312 FORMAT ( 30H IMPELLER BLADE EXIT BACKSwEEP *5X ♦ 23HI MPELLER BLADE EXIT 
1 HT • . 5X.32HVANELESS DIFFUSER DIAMETER RATIO, 5X ,28HVANE0 DIFFUSER S 
PETTING ANGLE) 

313 FORMAT ! 33X *33H TOTAL VANED DIFFUSER THROAT AREA, 5X .25HVANED OIFFUS 
1ER AREA RATIO) 

314 FCRMAT i 47X* 36H INLET CRITICAL VELOCITY RATIO ARRAY) 

315 FORMAT! 7X . F 7. 0 . 2 1 X.F6. 2. 23X .F 8. 7 ,24X , F 3. 1* 16X.F6.2) 

316 FORMAT! 7X . F3. 0, 20X.F4. 3. 1 5X , I 1 ,9X, F6. 0 , 1 7X.F4. 2 ,26X, F4.2 ) 

317 FORMAT 124X.F4. 1 .2 7X.F6.4.34X.F6.4) 

318 FORMAT! 13X.F4. 1 .26X.F6.5,27X,F6.3*30X.F5.2) 

319 FCRMATI47X.F7.6.28X.F4.2) 

320 FORMAT ( 8F 10.2 ) 

400 FORMAT! IX, 24H I RR AT IONAL EXIT TRIANGLE) 

40 2 FORMAT! IX .2 1FV0VCR ARRAY TOO LARGE) 

403 FCRMAT 1 IX . 34HWT FLOW LESS THAN SURGE FOR V/VCR=»F4 .2) 

FNC 
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Subroutine LININT(X1 , Y1 , X , Y , TN, MX, MY , F) 


This subroutine interpolates a value of maximum pressure recovery coefficient C* * 

It 

from a table of vaned diffuser throat Mach numbers and aerodynamic blockage given as 
input. 

XI input, vaned diffuser throat Mach number 
Y1 input, throat aerodynamic blockage 

X input array of throat Mach numbers 
Y input array of throat blockages 

TN input two-dimensional array of C* * corresponding to throat Mach numbers and 
blockages 

MX input, number of throat Mach numbers 
MY input, number of throat blockages 
F output, interpolated value of C* * 


SIBFTC L IN INT 

SUPROWT INE L IN (NT ( X l » VI .X , Y * TN. NX »MY , F I 
D IMFNS ICN X I MX > « V(MY>, TNI MX t MY I 
DC 1 .) 3=2. MX 

1 IF (XI .LE.XI J3) ) GO TO 2 

J3=NX 

2 00 3 J4=2,MY 

3 IF I YI .LE.YI J4I ) GO TO 4 

J4=MY 

4 Jl = J3- l 
J2=J4-1 

EPS 1 = 1 XI- XI Jl» W(XU3)-X(J11) 

FPS2=< Yl-YI J2» ) /( Y(J4>-Y( J2>> 

EPS3* 1 .-EPS1 
EPS 4= 1 .-FPS2 

F=TN( .1 1.J2 »*EPS3*EPS4+TNI J3, J2 >*EPS1*EPS4+TN 4 J 1, J4 >*EPS2*EPS3*TN( J 
13.J4)*EPS1*EPS2 
RETURN 
ENC 


Subroutine FNTGRL(NO, DELTAR, F, S) 


This subprogram is used to integrate the function 


M p (pa/p'a') in the vaneless dif- 


fuser . 
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NO number of equally spaced radii 

DELTAR radius ratio between stations 
F function M^(pa/p'a') 

S integral of F 

No FORTRAN listing of the subprogram is available but a brief description of the 
method of integration is as follows: 

The subprogram uses Simpson’s rule in the following manner: 

First interval: 



f(x)dx = 0 


Second interval: 



+ 8f 2 -f g )/12 


Remaining intervals: 



- f i-l 


)/12 


Lewis Research Center, 

National Aeronautics and Space Administration, 
and 

U.S. Army Air Mobility R&D Laboratory, 
Cleveland, Ohio, August 2, 1973, 

501-24. 
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APPENDIX A 


SYMBOLS 

A 2 

A area, m 

AR area ratio 

a local acoustic velocity, m/sec 

B effective depth ratio 

b blade height or passage depth, m 

C wall curvature, m~* 

C f skin friction coefficient 

C specific heat, J/(kg)(K) 

r 

C** maximum pressure recovery coefficient at a constant area ratio 
P 

D diameter, m 

Df diffusion factor 

e g inlet guide vane loss coefficient 

h spacing for numerical integration 

Ah incremental compressor specific enthalpy, J/kg 

K constant 

L blade length, m 

M Mach number 

N rotative speed, rpm 

n distance along radius, m 

o 

p pressure, N/m 

Q volumetric flow rate, m /sec 

q dimensionless enthalpy 

R gas constant, J/(kg)(K) 

R radius ratio 

Re Reynolds number based on inlet total conditions 

r radius, m 

s arc length, m 
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As incremental arc length, m 

T temperature , K 

t blade thickness, m 

u blade speed, m/sec 

V absolute gas velocity, m/sec 

W relative gas velocity , m/sec 

w mass flow rate, kg/sec 

Z number of impeller blades 

a absolute flow angle, deg from meridional plane 

relative flow angle, deg from meridional plane 
blade angle, deg from meridional plane 
y specific heat ratio 

6 ratio of inlet total pressure to standard sea-level pressure 

A5* incremental boundary layer displacement thickness, m 

e difference between compressor inlet relative flow angle and optimum incidence 

angle 

£ vaneless diffuser loss coefficient 

■q total efficiency 

Arj total efficiency decrement 

8 ratio of inlet total temperature to standard sea-level temperature 

A. inducer hub -tip diameter ratio 

o 

p dynamic viscosity, N-sec/m 

o 

v kinematic viscosity, m /sec 

p gas density, kg/m 

<p flow coefficient 

X wall slope or streamline angle, deg from axial 

Subscripts: 

act actual 

AD adiabatic 

aero aerodynamic 
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av 

BL 

calc 

choke 

cr 

DES 

DF 

eff 

eq 

est 

exit 

H 

HYD 

i 

id 

IGV 

INC 

is 

L 

m 

max 

MF 

n 

opt 

R 

RC 

SF 

SL 

st 

surge 


average 
blade loading 
calculated 

value at choking flow rate 

critical condition 

design 

disk friction 

effective 

equivalent 

estimated 

exit 

hub 

hydraulic 

iteration subscript 
ideal 

inlet guide vane 
incidence 
isentropic 
lost 

meridional 

maximum 

rms 

normal 

optimum 

rotor 

recirculation 
skin friction 
slip 

stagger 

value at surge flow rate 
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T tip 

u tangential 

v vortex 

VD vaned diffuser 

VLD vaneless diffuser 

0 station just upstream of inlet guide vanes 

1 impeller inlet upstream of blade 

2 impeller inlet just downstream of blade leading edge 

3 impeller exit 

4 vaned diffuser leading edge 

5 vaned diffuser throat 

6 vaned diffuser exit 
Superscripts : 

’ absolute total condition 

” relative total condition 
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APPENDIX B 


EQUATIONS 

This section presents a more detailed discussion of the method of calculating com- 
pressor performance than was presented in the section Method of Analysis. 

The equations used in calculating the flow solution through the compressor are listed 
here in the approximate order of solution in the FORTRAN pregram. The known con- 
stants are the input fluid state and thermodynamic properties, compressor geometry, 
and velocity diagram data. The working fluid properties and state conditions are spe- 
cific heat ratio y , inlet total pressure p^, inlet total temperature T|j, and inlet dy- 
namic viscosity based on total conditions The velocity diagram data are compressor 

design rotative speed N^gg, the fraction of design rotative speed for which performance 
data is desired N/N DEg , the number of values of (V/V cr ^ for which performance is to 

be calculated, the values of (v /V cr ) to be used, and the inducer solid -body swirl angle 

at the rms diameter The compressor data are inducer inlet hub and tip wall 

slopes Xjjj Xj-p? inducer inlet hub and tip wall curvatures C^g and C^rp, inducer 

inlet tip diameter D^, inducer inlet hub-tip diameter ratio X, inducer -inlet-tip to 
impeller -exit diameter ratio D^/D^, number of blades at impeller exit (including split- 
ters) Zg, impeller blade exit angle (positive opposite to direction of rotation ) im- 
peller exit blade height bg, inducer inlet blade angle at the rms diameter in- 

ducer inlet blockage factor calculated at the rms diameter, vaneless diffuser diameter 
ratio D^/Dg, vaned diffuser setting angle a 4, total vaned diffuser throat area A^, 
vaned diffuser area ratio Ag/A*., and wetted surface friction coefficient C^. With these 
known values compressor performance can be evaluated by the method described. 


Calculation of Inlet Velocity Triangles and Compressor Weight Flow 

Calculation of the inducer inlet velocity diagrams and weight flow are handled by two 
methods depending on whether the inlet is swirl-free or subject to solid -body swirl. For 
the swirl cases it is assumed that the inlet guide vanes are located in a constant area 
annulus with zero wall slope with respect to the axial direction. The inner and outer 
diameters of the annulus are assumed to be equal to the inducer inlet hub and tip diam- 
eters. For the swirl-free cases it is assumed that the flow orthogonal can be approxi- 
mated by a radial line. These assumptions greatly simplify calculation of the velocities 
specified by the general velocity gradient equation. 
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Swirl-free inlets . - For compressor configurations with swirl-free inlets the input 
value of /V cr j is"used as the hub streamline value of V/V cr - The variation of inlet 

velocity with radius is calculated from the general velocity gradient equation simplified 
for swirl-free flow. This equation is (in simplified form from ref. 10) 


V m = V mlH 



(Bl) 


where V m ^jj is the hub absolute velocity, C is the streamline curvature, and n is the 
distance along a radius. Streamline curvature and slope at the rms inlet diameter are 
approximated by linear interpolation of the input hub and tip values. The integral in the 
previous equation is approximated by the trapezoidal rule at the rms diameter and by 
Simpson’s rule for unequal intervals at the tip. The numerical approximations are 



ih +c imf) 


V mlMF " V mlH e 


(B2) 


and 


' h o +h D 


V mlT “ V mlH e 



( V h i) 

c + u c + 

V 1H h i h o 1MF 



(B3) 


where 


h 0 ~ r lMF _ r lH 


(B4) 


and 


h l ' r lT " r lMF ( B5 ) 

Static densities at the hub, rms, and tip diameters are determined from the inlet total 
density and the calculated velocity at each station. Velocity components normal to the 
annular inlet plane are calculated using the velocities and streamline angles at the three 
stations: 
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V mln = V ml c ° s X < B6 » 

where x is the streamline angle. Inlet mass flow is then computed using Simpson's 
rule for unequal intervals 


w = 2tt 


h„ + h. 


“ ( 2 Vp r) lH V mlHn + 


( h o + h i) : 


IMF v mlMFn 



Inlets with solid -body swirl . - For impeller inlets with inlet guide vanes producing 
solid-body swirl the equations presented in reference 2 for kinetic energy loss were 
iteratively solved for state conditions and velocities downstream of the inlet guide vanes 
Conditions at the rms diameter were assumed to be representative of the flow distribu- 
tion at the impeller inlet. The inlet guide vane stagger angle was assumed to be one- 
half of the turning angle at the rms diameter; that is, 



0.0076 

cos " 0.025 


COS Oft 


Using the input value of (v/V cr j as a first approximation to the value of V/V cr down- 
stream of the inlet guide vanes at the rms diameter, the following equations are solved 


iteratively with successive increments in V/V cr until continuity checks the calculation 
of equation (B7) : 
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2 


cr 


KE 


cr 


'A - n - 

VP '/IMF 


KE 
C T' 


.y/(y-l) 


KE id = 


KE 

1 - e c 


P 1MF , / 1 .“ld Xr/6 '' 1) 


Po 


Vi), 


P 1MF 


Pq 

P 1MF “ p 0 


p Amf 


P 1MF 


P 1MF 


RT’ 


P 1MF " P 1MF ^ j 


IMF 


P 1MF ' P 1MF fe 


1/V 


IMF 


Q 1 “ ^ D 1T^ “ (~~ ) V cr cos “IMF 


cr. 


IMF 


w = Q l p lMF 


(BIO) 

(Bll) 

(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 

(B18) 

(B19) 
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When continuity is satisfied downstream of the inlet guide vanes, the impeller inlet 
velocity triangles are calculated. Combination of the Euler equation with the solid-body 
vortex equation and simple radial equilibrium equation shows that the meridional veloc- 
ity distribution can be expressed as 


V 


ml “ 



(B20) 


Using the value of V/V cr which satisfies continuity downstream of the inlet guide vanes 
and the input value of the rms turning angle allows solution of the constant in equa- 
tion (B20) . The tangential components of the inlet absolute velocity at the hub and tip are 
calculated from the vortex condition. Then the hub and tip meridional velocities can be 
determined as follows: 


V u = K Vu r 


(B21) 


where 


K 


Vu 


' ulMF 
r lMF 


(B22) 


The following general inducer inlet velocity diagram characteristics are calculated 
after the particular swirl-free or solid -body swirl calculations are completed: 


V 


IT 




y2 + y2 

v mlT + v ulT 


(B23) 


— sin 


• -1/ V ulT 


IT, 


(B24) 


W ulT “ U 1T " V ulT 


(B25) 


0it 


= tan 


-1/WulT 

r mlT/ 


(B26) 
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W 1T “ V V mlT + W ulT 


W ulH “ U 1H ' V ulH 


*>1H 


= tan 


-1/ W ulH 


mlH, 


W 1H= / V mlH +W ulH 


IMF 




2 2 
ulMF + V mlMF 


W ulMF ” U 1MF “ V ulMF 


W 


IMF 




v mlMF ulMF 


%MF 


= tan 


-If W ulMF 
^mlMF/ 


(B27) 

(B28) 

(B29) 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 


The wheel speeds used in the previous equations are obtained from the following input 
constants: 


rr | — 1 N DES D 1T 


N. 


u 


DES/ 


IT 


60 


U 1H " XU 1T 


D 


U 1MF “ U 1T 


IMF 


D 


IT 


(B35) 

(B36) 

(B37) 
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where 


D 


IMF 


vF 


j T (l + A 2 ) 


(B38) 


Individual Losses 


Inlet guide vane loss . - When the inducer inlet velocity diagram characteristics 
have been determined by the method described in the section Inlets With Solid-Body 
Swirl, the inlet guide vane loss is computed using the equation found in reference 1, sub- 
ject to the assumptions made in reference 3. 


Ah IGV " e S KE id 


(B39) 


Inducer incidence loss . - The equations developed in reference 4 are used to deter- 
mine the optimum inducer incidence angle for incompressible flow. Compressibility 
effects are ignored because the incidence loss is evaluated at the rms diameter and com- 
pressibility effects are insignificant except at the hub according to reference 4. Using 
the input blade blockage factor at the rms diameter the optimum incidence angle is found 
by 


€ = tan 


-1 


(1 - B 1 )tan g 1MF 
1 + Bj tan 2 0 1MF 


(B40) 


where 


ZA 

Bj = 1 (B41) 

27rr lMF 


The optimum relative flow angle is then 


^opt ' %MF “ 6 < B42) 


and the component of the relative velocity lost is 
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W L " W 1MF sin l^opt ' ^1MF 


The enthalpy loss due to incidence is expressed as 


w; 


Ah 


INC 


2C, 


and the total pressure just inside the blade row is 


P 2MF “ p 'lMF 


e 




(B43) 


(B44) 


(B45) 


Impeller work and losses . - The impeller exit density is obtained through iteration 
of the impeller loss equations and calculated state properties. The enthalpy rise in the 
impeller is initially approximated as 


Ah est = u 3 


(B46) 


and the corresponding impeller exit total temperature is then 


'Ah 


est 


36St l p rpt / 0 

So 


(B47) 


The impeller exit density is initially approximated to be 


/ T wV /fy-1) 

p 3 " PimfI “^T 


(B48) 


The meridional component of exit absolute velocity is calculated from continuity 


w 


m3 


7rp 3 D 3 b 3 


(B49) 
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The tangential component of exit absolute velocity is then 


V u3 = u 3 - V m3 tan %b - V SL 


where the slip velocity Vg^ is calculated by the method of reference 11. 

_ u 3 V C0S ?3b 


SL 


,0.7 


The remaining state and flow properties at the exit were calculated from the 
ing relations: 


npM _ rp, t I 
1 3 “ IMF + 


u 3 ~ U 1MF 
2C„ 


W u3 = u 3 - V u3 


W 3= V^mS^uS 


w; 

rp __ rp t 1 O 

3 3 "^c" 


V 3 = y V L3 + V u3 


T 3 = T 3 + 


2C 


The work loss due to blade loading is calculated from the equation 


ih BL = 0.05D f 2 u| 


which is given by reference 5. The diffusion factor Dj is defined as 


(B50) 

(B51) 

follow- 

(B52) 

(B53) 

(B54) 

(B55) 

(B56) 

(B57) 

(B58) 
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W 3 + K BL q aero 



- 

- 

W 1T 

Z 3 L _ d 1t'i 

^1T 
+ 2 — 

u 3 

1 

o 

CO 

1 

CO 

Q 


(B59) 


where a value of 0.75 is used for Kg^ for conventional impellers and a value of 0.6 is 
used for impellers with splitters. A parametric study of calculated diffusion factors 
with a variation in the number of blades indicated that changing the constant to 0.6 would 
compensate for the changing solidity near the exit. The dimensionless actual head is 
obtained by dividing the enthalpy rise by the exit blade speed squared 


q 


aero 


Ah 


aero 


u 


2 

3 


where 


Ah 


aero 




(B60) 


(B61) 


Impeller disk friction loss is computed using the method of reference 6 


Ah DF = 0.01356 


P 3 


wRe 


0.2 


" 3 D 3 


(B62) 


The Reynolds number in this equation is based on impeller exit dimensions and inlet 
total conditions 


Re 


u 3 D 3 p lMF 

p'o 


(B63) 


Skin friction loss is calculated from the correlation found in reference 5 
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Some studies of the effects of the addition of splitter blades on the impeller flow in- 
dicated that the mean channel relative velocity is increased in the splitter region. The 
ratio of mean flow path blade length and hydraulic diameter -exit diameter ratio remain 
virtually unchanged. Computation of the mean relative velocity ratios for an impeller 
both with and without splitters showed that, for splitters, the average relative velocity is 
about 14 percent higher than for conventional impellers. This phenomenon is accounted 
for by using the constant 7.0 for impellers with splitters and 5.6 for impellers without 
splitters. 

At this point in the iteration procedure the impeller ideal enthalpy rise is calculated 
as follows: 
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Ah id = Ah aero “ Ah IGV “ Ah INC “ Ah SF " Ah BL 


(B68) 


Defining the impeller efficiency as 




Ah id 


Ah. 


aero 


the state conditions at the impeller exit are calculated from the relations 


Ah aero 

P3-P2MFI . T , 

l p 0 


,r /(y -l) 


p, = 


Pi 




^3 



(B69) 


(B70) 


(B71) 


(B72) 


The value of pg obtained from equation (B72) is compared to the value estimated in 
equation (B48) . 

If the values do not compare within 0. 1 percent, the former value is used as a new 
approximation and equations (B49) to (B72) are iterated until two consecutive values of 
exit density agree within the tolerance specified. 

Recirculation loss . - Recirculation of the working fluid from the vaneless space back 
into the impeller results in additional work which is lost to the system. Reference 5 ex- 
presses this loss as a function of the impeller exit absolute flow angle, impeller diffu- 
sion factor, and impeller exit tip speed. The following modified form of the equation of 
reference 5 is used to calculate the recirculation loss: 

Ah RC = 0,02 y <*3 D f«3 (B73) 

Vaneless diffuser loss . - Reference 7 developed differential equations relating Mach 
number, flow angle, and total temperature to vaneless diffuser radius ratio through the 
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fundamental relations of continuity, equilibrium, heat transfer, and fluid state. These 
equations are evaluated at ten equally spaced intervals in the vaneless space. Two iter- 
ations of the equations are performed at each station; the first iteration uses the calcu- 
lated properties at the point of interest to compute the variation in Mach number and 
flow angle to the next point; the second iteration uses the estimated properties of the 
downstream point to compute the same variations. The arithmetic average of the two 
iterations is then used for the solution. 

These equations, simplified for adiabatic flow in a geometrically constant depth 
radial passage are 


-2(1 +21^1 m 2 ' ■- 


1 dM* 


M 2 dR 


M 2 - sec 2 a 


(yM 2 - tan 2 a) 


1_ dB _ sec a 


B cos of B 


dR R 


(B74) 


and 


1 d tan a 


sec 


ten a dR m 2 - sec 


— — /[l + (y - 1)M 2 1 -S _ + l<g-g 

;2 a -I B cos a B dR R 


A 


(B75) 


where 


? = 


C f r 3 


(B76) 


R = — 
r 3 

B = — 


(B77) 


(B78) 


The schedule of effective passage depth is determined by boundary layer displace- 
ment thickness growth on the end walls. The end walls in the vaneless diffuser are as- 
sumed to be parallel and have a spacing equal to the impeller exit blade height. The 
flow pattern between the end walls is approximated as a log spiral between adjacent cal- 
culation stations. With these assumptions the incremental flow path length is 
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(B79) 


AS = 


r 2 AR 
COS O' 


and the incremental boundary layer displacement thickness is 


/V \“®* ^ 

A8* = 0.037 S" 0 - 2 — I AS 


where 

n 

s - E AS > 

i=l 

The effective passage depth is 


B i + 1 = B i * 


2 A6* 


and 


AB = B. - B. 1 
l l+l 


(B80) 


(B81) 


(B82) 


(B83) 


Loss the total pressure in the vaneless space is computed from the following equa- 
tion derived in reference 4: 


1 




COS Qfg 



(B84) 


The integral in equation (B84) is determined by numerical methods. Vaned diffuser 
leading edge static pressure is calculated from the isentropic relation using the total 
pressure obtained from equation (B84) and the Mach number obtained from the distribu- 
tion prescribed by equation (B74) 
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(B85) 


P 4 = 


p 4 


1 +- 




Vaneless diffuser loss is then determined by 


Ah VLD " C p T 3 1 



(r-l )/V 



(y-l)/Vl 


(B86) 


Vaned diffuser loss . - Vaned diffuser performance is predicted by use of the test 
data reported in reference 8. Lines of maximum pressure recovery coefficient at a 
given area ratio were estimated from the performance maps reported for single plane 
divergence diffusers with square throats. The component of the vaned diffuser leading 
edge Mach number parallel to the vane setting angle is used as the throat Mach number. 

A loss in total pressure associated with this incidence is calculated holding static pres- 
sure constant between the leading edge and throat. The throat blockage is estimated 
from the displacement thickness growth on the vaneless diffuser end walls. Using these 
conditions of throat Mach number, total pressure, and blockage the pressure recovery 
of the vaned diffuser is extrapolated for the specified geometric area ratio. Vaned dif- 
fuser exit static pressure is then calculated from the relation 

P6 = C p* (p 5-P5» + P5 (B87) 

The exit total pressure and velocity are determined by trial and error using esti- 
mated values of exit critical velocity ratio. The total pressure at the vaned diffuser 
throat is assumed to be constant across the free-stream area outside the boundary layer 
displacement thickness. The total pressure and velocity at the vaned diffuser exit are 
calculated assuming full flow across the geometric area. Starting with a low value of 
( V /Vcr) l° wer estimate of exit total pressure is expressed as 



p 6 


■1-TLZ± 

y + 1 



(B88) 
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and the higher estimate is 


where 




(B89) 


(B90) 


Equations (B88) and 


(B89) are solved using increasing values of 



until conver- 


gence is achieved. The vaned diffuser loss is then calculated from the relation 


Ah VD " C p T 3 



(y-1 )/y 



(y-D/y I 


(B91) 


Efficiency decrements . - Decrements in efficiency caused by the individual losses 
are obtained from the following equations: 


Hgv 


At? inc 


At >bl 


Ah 


IGV 


Ah 


act 


Ah 


INC 


Ah 


act 


Ah 


BL 


Ah 


act 


(B92) 


(B93) 


(B94) 


At ? S F ~ 


Ah 


SF 


Ah 


act 


(B95) 
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(B96) 


Ah 


DF 


Ah 

Ah 


DF 

act 


i7 hic 


Ah. 

Ah 


RC 

act 


Atj 


Ah 


VLD 


VLD 


Ah 


act 


At;. 


VD 


Ah 

Ah 


VD 

act 


where 


Ah act “ Ah aero + Ah RC + Ah DF 


Overall efficiency . - The overall total efficiency is calculated as follows: 


^AD 


_ ^aero “ ^ Ah IGV + Ah INC + Ah BL + Ah SF + Ah VLD + Ah VD^ 


Ah aero + Ah RC + Ah DF 


(B97) 


(B98) 


(B99) 


(B100) 


(B101) 
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